Measuring the Out-of-Equilibrium Splitting of the Kondo Resonance 
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An experiment is proposed to measure the out-of-equilibrium splitting of the Kondo resonance in an 
ultrasmall quantum dot, by adding a third, weakly coupled lead to the standard two-lead quantum- 
dot system, and sweeping the chemical potential of that lead. Fixing the voltage bias between the 
source and drain leads, we show that the differential conductance for the current through the third 
lead traces the out-of-equilibrium dot density of states (DOS) for the two-lead system. This enables 
one to measure the dot DOS in the presence of an applied voltage bias. We show that this method 
is robust, and extends also to the case where the coupling to the third lead is no longer weak. 



PACS numbers: 72.15.Qm, 73.23.Hk, 75.20.Hr 

The understanding of strong electronic correlations 
far from thermal equilibrium is one of the challenging 
problems in contemporary mesoscopic physics. A pri- 
mary example is the out-of-equilibrium Kondo effect, re- 
cently measured in ultrasmall quantum dots . In the 
Kondo effect, an impurity moment undergoes a many- 
body screening by the conduction electrons, producing a 
narrow low-energy resonance in the impurity density of 
states (DOS) at low temperatures ||. This Kondo reso- 
nance is responsible for the well-known enhancement of 
the low-temperature resistivity in dilute magnetic alloys, 
and for the enhancement of the low-temperature conduc- 
tance in tunneling through small quantum dots with an 
odd number of electrons. One of the striking predictions 
for the out-of-equilibrium Kondo effect in quantum dots 
is the splitting of the Kondo resonance in the dot DOS 
for a finite applied bias . 

Although well-established theoretically, it remains un- 
clear whether one can measure the above splitting in the 
dot DOS. The splitting of the Kondo resonance does not 
show up in the differential conductance for the current 
through the dot since, unlike in conventional noninter- 
acting systems, the DOS itself is strongly voltage depen- 
dent. Photoemission, besides being limited in resolution, 
measures the "occupied" DOS, i.e., the product of the 
DOS and the effective distribution function. Since the 
latter distribution is unknown away from thermal equi- 
librium, one can not reliably extract the underlying DOS 
using photoemission. A key assumption in conventional 
tunneling measurements of the DOS is that the weakly 
coupled sample and probe are each effectively in equilib- 
rium. Clearly this assumption breaks down for the quan- 
tum dot, which raises the fundamental question: Can one 
actually measure the DOS of an interacting mesoscopic 
system far from thermal equilibrium? 

In this paper we show that the out-of-equilibrium split- 
ting of the Kondo resonance can indeed be measured, by 
adding a third, weakly coupled lead to the standard two- 
lead quantum-dot system jffl. Fixing the voltage bias 




FIG. 1. A schematic sketch of the proposed apparatus. 
An ultrasmall quantum dot is coupled by tunneling to three 
metallic leads, each of which is kept at a separate chemical 
potential. The corresponding tunneling matrix elements are 
controlled by varying the potential barriers, while the dot 
energy level is adjusted by applying a gate voltage. The dot 
level and the couplings are to be tuned such that the dot is in 
the Kondo regime, and the coupling to the third lead is much 
weaker than to the other two leads. Fixing the source-drain 
voltage bias at [4% — \i\ — eV s d and sweeping the chemical 
potential [13 by varying V31, one measures the current ^(Viai) 
between the dot and the third lead. Up to thermal broadening 
and rescaling, the differential conductance G3 = dlz/dVsi 
traces the out-of-equilibrium two-lead dot DOS. 



between the source and drain leads and sweeping the 
chemical potential of the third lead, we show that the 
differential conductance for the current through that lead 
traces the out-of-equilibrium two-lead dot DOS, much in 
the same way as in conventional tunneling measurements 
of the equilibrium DOS. In this manner, one can measure 
the dot DOS in the presence of an applied voltage bias. 
Furthermore, we find that the splitting of the Kondo res- 
onance can be observed also at temperatures above the 
Kondo temperature, and that the same qualitative fea- 
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tures are obtained when the coupling to the third lead 
is no longer weak. Contrary to a previous proposal to 
probe the splitting of the Kondo resonance using two ca- 
pacitively coupled quantum dots S , our approach allows 
for a direct measurement of the nonequilibrium DOS. 

Figure |l] shows the proposed experimental setting. It 
consists of an ultrasmall quantum dot, coupled to three 
metallic leads. Each lead is kept at a separate chemi- 
cal potential /j, a (a — 1, 2, 3), and is coupled to the dot 
via a tunneling matrix element t a . The energy level on 
the dot, Ed, is controlled by applying a gate voltage. 
Modeling the charging energy on the dot by a Hubbard 
on-site repulsion U, the system is described by a gener- 
alized three-lead Anderson impurity model: 



^=EE [( £fe + V a )c X aka C aka + t a (J aka da + /i.e.) 



a—1 



Here c„ u „ creates an electron with wave number k and 
spin projection a in lead a, while S a creates a dot electron 
with spin projection a. In the following we assume that 
the tunneling matrix elements and the energy level of the 
dot are adjusted in such a way that a stable local moment 
is formed on the dot. This corresponds to the condition 
r = ELi r " < -E D ,E D + U, where T a = p a irt 2 a is 
the hybridization width associated with lead a (p a is the 
conduction-electron DOS at the Fermi level of lead a). 
We shall generally regard T3 as much smaller than Ti and 
r 2 , though comparable couplings will also be considered. 

The main quantities of interest are the electrical cur- 
rents flowing from the dot to each of the three leads. 
The electrical current flowing into lead a is given by 
the nonequilibrium average of the current operator I a — 
it a (e/h) ^2 ka {c^ aka d a — h.c.} (— e is the electron charge). 
This average is conveniently expressed via the lesser and 
retarded Keldysh Green functions for the dot electrons: 



+E D ^2 d l d <y + Ud\d^d\di. 



(1) 



= / (dUAt))e ltt dt, 



(2) 



G da {e) = -z {{d a (t),dt})e^dt. (3) 
Jo 

Specifically, generalizing the expression of Meir and 
Wingreen || from two to three leads one obtains 

4e f°° 

(4) 

where A d (e) = — -^Im{G d (e)} is the dot spectral function 
(i.e., DOS); /(e) is the Fermi-Dirac distribution function; 
and u a (e) = p a (e)/p a is the reduced conduction-electron 
DOS in lead a. Here we use the notation by which the 
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FIG. 2. Comparison of the out-of-equilibrium dot DOS 
for a two-lead system (full lines) with the differential conduc- 
tance 03(^,3) in the proposed three-lead setting (empty cir- 
cles). All calculations were performed using the NCA. Here 
Ti = T 2 = D/30 (symmetric coupling), E d /D = -0.278, 

ryih = 0.01, t = T K , and Go = 4e 2 r 3 (r! + r 2 )/Dhr. 

The two-lead dot DOS is scaled by the bandwidth D. The 
source-drain voltage bias, /12 — /ii = eV s d, takes the values 
eVsd/ksTK = 0,10,25,50, and 100, as indicated by the at- 
tached labels. Up to thermal broadening, the scaled differ- 
ential conductance coincides with the out-of-equilibrium dot 
DOS, in accordance with Eq. (|^). Inset: An extended im- 
age of the dot DOS for a two-lead system with T = Tk and 
eV sd /k B T K = 50. 



arguments of p a (e) and v a (e) are measured relative to 
the corresponding chemical potential, /i Q . Note that we 
have restricted ourselves in Eq. (^) to a zero magnetic 
field, in which case all spin dependences drop. 

The integration over energy in Eq. (^) contains a nat- 
ural cut-off, namely, the maximal chemical-potential dif- 
ference between the leads, broadened by the tempera- 
ture. To the extent that one can neglect the energy de- 
pendence of v a (e) on that scale, it is possible to elim- 
inate G d from the expression for ^3 using the identity 
h = J 3 (Ti + r 2 )/r - (h + I 2 )T 3 /T: 



h = 



4er ; 



A d (e)[(T 1 + r 2 )/(e-/i 3 ) - ^/(e-Mi) 



-r 2 /(e-/i 2 )] de. 



(5) 



Analogous expressions apply to each of I\ and l-i. 

Equation (|^) generalizes the standard two-lead expres- 
sion for tunneling through an Anderson impurity It 
retains the familiar form of an integral of a spectral func- 
tion times a weighted difference of Fermi functions. If 

(e) were insensitive to variations in the chemical po- 
tentials (as is the case for conventional noninteracting 
systems) , then the derivative of ^3 with respect to any of 
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the Ha's would have given the thermally broadened dot 
DOS at fi a . Clearly, this is not the case for a quantum 
dot in the Kondo regime. Nevertheless, for T 3 <C Ti +T 2 
one can still fix the chemical potentials [i\ and [ii and dif- 
ferentiate the current I3 with respect to [13 , to obtain the 
out-of-equilibrium two-lead dot DOS for a source-drain 
voltage bias of [ii — [i\ = eV s d- To see this we note that, 
for T3 <C Ti + T2, the impurity Green functions only 
weakly depend on the coupling to the third lead. Hence 
A d (e) and G^(e) are insensitive to variations in /13, de- 
spite the strong dependence on \i\ and \i2- Fixing \i\ and 
[i2 and sweeping the bias V31 (see Fig. |l|), the differential 
conductance (^3(^3) = dl^/dV-n reduces then to 



4e^3 



(ri + r 2 ) 



A d (e+» 3 )^de. (6) 
oe 



Hence, up to rescaling, G^ip^) is just the thermally 
broadened spectral function ^(^3), which in effect is 
the two-lead dot DOS for the chemical potentials /ii and 
/i2- This means that the weakly coupled third lead acts 
as a conventional tunneling probe for the dot DOS, even 
though the dot itself is far from thermal equilibrium. 

To test the above result, we have computed the dif- 
ferential conductance 6*3(^3) by direct numerical dif- 
ferentiation of the current I3 of Eq. (|J) with respect 
to V31, without resorting to Eqs. (||) and (||). Focus- 
ing for convenience on the limit U — ► 00 (double occu- 
pancy is forbidden on the dot), we evaluated the impurity 
Keldysh Green functions using the noncrossing approx- 
imation (NCA) ffl. The NCA is a self-consistent per- 
turbation theory about the atomic limit, which is known 
to provide a good quantitative description of the tem- 
perature range T > Tk- It has been extensively used 
to study dilute magnetic alloys and has been suc- 
cessfully applied to the out-of-equilibrium Kondo effect 
for both single-channel ||[ll],[lj] and two-channel scat- 
terers [O. In the present context it has the crucial ad- 
vantage that it can be easily generalized to a multi-lead 
setting. As a large- N theory, though, the NCA fails to re- 
cover Fermi- liquid behavior at low temperatures ||l0|| , and 
it overshoots the unitary limit for T < Tk in the N = 2, 
nondegenerate case. To avoid these NCA pathologies, we 
restrict attention hereafter to the range T > Tk- 

Figure || compares the the out-of-equilibrium dot DOS 
for a symmetrically coupled two-lead system with the 
differential conductance 6*3(^3) in the proposed three- 
lead setting. Here and throughout the paper we set 
(/ii + /^2)/2 as our reference energy, and use a semi- 
circular conduction-electron density of states with half- 
width D: v a {e) = y/l - (e/D) 2 . The impurity model 
parameters are E d /D = -0.278 and T x = T 2 = D/30, 
corresponding to a Kondo temperature of k B TK / D = 
2.5xl0~ 4 @. The third-lead coupling, r 3 /Fi = 0.01, is 
sufficiently weak as not to affect Tk- 

Up to thermal broadening, the scaled differential con- 
ductance of Fig. coincides with the out-of-equilibrium 
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FIG. 3. The out-of-equilibrium two-lead dot DOS (full 
lines) versus the three- lead differential conductance 03(113) 
(empty circles), for a fixed source-drain voltage bias of 
eVsd/ksTK = 50, and different ratios of Ti to IV Here 
(ri + T2)/D = 0.067 is kept fixed as to maintain the same 
Kondo temperature, ksTx/D — 2.5 x 10~ 4 . The remain- 
ing parameters are E d /D = -0.278, F 3 /(ri + Y 2 ) = 0.005, 
T = T K , and Go = 4e 2 r 3 (ri + T 2 )/DhT. The ratio Ti/r 2 
takes the values 1/3, 1, and 3, according to the labels attached. 
The effect of an asymmetry in Fi and T2 is to enhance (re- 
duce) the peak at the chemical potential of the more strongly 
(weakly) coupled lead. Inset: Temperature dependence of 
03(^3), for eVsd/ksTK = 50. Full, dashed and dotted lines 
correspond to T /Tk = 1,5, and 10, respectively. The effect 
of a temperature is to smear the peak structure of 03(^3). 



two-lead dot DOS, in excellent agreement with Eq. (g). 
For zero source-drain voltage bias, eV s d = 0, there is 
a single Kondo peak in G 3 {^jlj,), corresponding to the 
thermally broadened Abrikosov-Suhl resonance in Ad(e). 
[Note that the height of the resonance in Ad(e) is some- 
what lower than 1/irT, indicating that the Kondo effect 
is not yet fully developed.] As the source-drain bias is 
increased, a two-peak structure gradually develops in 
G 3(113). Two well-resolved peaks are finally observed 
once eV s d exceeds IzbTk and ksT by an order of mag- 
nitude or more. For eV s d 3> fcsT, fc^T^-, the effect of 
thermal broadening is quite negligible, as the widths of 
the two Kondo peaks are governed by the dissipative life- 
time induced by V sd ||, rather than by ksT or ksTK- 
The asymmetric Kondo line shapes, particularly for small 
to intermediate values of eV s d, stem from the absence of 
particle-hole symmetry in the infinite-?/ Anderson model. 

Figure || depicts the out-of-equilibrium two-lead dot 
DOS versus G 3 (/x 3 ), for eV sd /k B T K = 50, T = T K , and 
different ratios of r x to r 2 . The couplings (Fi + T 2 ) 
and r 3 / (ri + T2) are kept fixed as to maintain the same 
Tk- As before, excellent agreement is obtained between 
^3(^3) and the two-lead dot DOS, in accordance with 
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FIG. 4. The three-lead difFerential conductance Gs(/i3), 
for ri = r 2 and different ratios of F3 to Ti- The source-drain 
voltage bias, eV s d/kBTK, is equal to zero for the three 
single-peak curves, and 50 for the three split-peak curves. 
T/D = 0.067 is kept fixed in all curves as to maintain the 
same Kondo temperature. Here Ed/D — —0.278, T = Tk, 
and Go = Dftr/8e 2 r 3 (ri + T 2 ). The overall peak structure 
of G 3(^13) is mostly unchanged when T3/T is increased. Right 
inset: The evolution of Gs(/i3) with increasing source-drain 
voltage bias, for I\ = T 2 = T 3 and T = T K - Here eV sd /k B T K 
equals 0,25,50, and 100. Left inset: The dot DOS, DA d (e), 
for eV sd /k B T K = 50, fi 3 /k B T K = 0, and I?i = T 2 = T 3 . 
There are three separate peaks in the dot DOS, one peak at 
the chemical potential of each of the three leads. 



Eq. (0). The effect of an asymmetry in the couplings to 
the source and drain leads is to enhance the peak at the 
chemical potential of the more strongly coupled lead, at 
the expense of the peak at the chemical potential of the 
more weakly coupled lead. As seen in Fig. |], the effect 
is quite dramatic. Only a shallow peak is left at the 
chemical potential of the more weakly coupled lead when 
the ratio between the two couplings is three. Accordingly, 
the dot DOS increasingly resembles that for sole coupling 
to just a single lead. The inequivalent line shapes for 
T1/T2 = 3 and r^/Ti = 3 are again the result of the 
absence of particle-hole symmetry in our model. 

The effect of a temperature is to smear the peak struc- 
ture of Gs(fi3), as seen in the inset of Fig. ||. This smear- 
ing stems both from the thermally broadened Kondo 
peaks in the dot DOS, and from the convolution with 
the derivative of the Fermi-Dirac function in Eq. (||). 
The two-peak structure of the Kondo resonance remains 
visible, though, up to UbT about an order of magnitude 
smaller than eV s d (see also Fig. ||). For large source-drain 
bias, as in Fig. |^, the two-peak structure thus extends up 
to temperatures well above Tk- 

Thus far our discussion has been restricted to weak 
coupling between the dot and third lead. One may ask 



how does the qualitative picture change once this cou- 
pling is no longer weak, especially as the reasoning lead- 
ing to Eq. (H) no longer holds. In Fig. [| we have plotted 
G^(n'i) for Y\ = T2 and different ratios of T% to Y\. Here 
r is kept fixed as to maintain the same Tk- Despite the 
notable changes in the dot DOS, GzijJ-z) remains mostly 
unchanged even for T 3 as large as Ti. Indeed, while A^(e) 
acquires a three-peak structure for ri = T2 = T3 (see 
Fig. [|, left inset), (^3(^3) remains closely related to the 
two-lead dot DOS for the chemical potentials [i\ and fi2- 
This clearly shows the robustness of the proposed method 
for measuring the splitting of the Kondo resonance. 
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olated zero-temperature limit. This definition of Tk is 
common in dilute magnetic alloys, where R(T) corre- 
sponds to the electrical resistivity due to scattering off 
the magnetic impurities. With Tk so defined, the half- 
width of the Abrikosov-Suhl resonance in the NCA im- 
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